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Abstract. The scattering cross section for a long- wavelength planar gravitational wave 
impinging upon a rotating black hole is calculated, for the special case in which the direction 
of incidence is aligned with the rotation axis. We show that black hole rotation leads to a term 
in the cross section that is proportional to aw. Hence, contrary to some claims, co-rotating 
and counter-rotating hclicities are scattered differently, and a partial polarization is induced 
in an unpolarized incident wave. 

The scattering amplitudes are found via partial wave series. To compute the series, two 
ingredients are required: phase shifts and spin-weighted spheroidal harmonics. We show 
that the phase shifts may be found from low-frequency solutions of the radial Teukolsky 
equation derived by Mano, Suzuki and Takasugi. The spheroidal harmonics may be expanded 
in spherical harmonics; we present expansions accurate to second order in auj. The two 
ingredients are combined to give explicit expressions for the helicity-conserving and helicity- 
reversing amplitudes, valid in the long-wavelength limit. 
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1. Introduction 



Gravitational waves (GWs) are propagating ripples in space-time whose existence is predicted 
by General Relativity. There is strong indirect evidence for their existence, for example, from 
thirty-five years of pulsar timing measurements [1] . Yet due to the tiny expected amplitude 
of waves reaching Earth, GWs have not been measured directly. Now, nine decades after 
the formulation of Einstein's theory, many experimentalists are optimistic that "first light" 
detections are imminent, at either (existing) ground-based [2] or (future) space-based [3] 
interferometers. 

GWs are of interest to astronomers because they are generated by some of the most 
energetic astrophysical processes, such as binary mergers, supernovae and galaxy collisions. 
Electromagnetic radiation carries relatively little information about the most energetic regions 
of such processes, because photons are strongly scattered, absorbed and thermalized by 
intervening matter. On the other hand, gravitational waves are only weakly coupled to 
matter, and carry information about the dynamics at the heart of such processes. 

In this note, we are motivated by a simple general question: does the rotation of matter 
induce a polarization in a gravitational wave? To explore this question, we study a special case. 
We compute the differential scattering cross section for a long-wavelength planar gravitational 
wave impinging along the rotation axis of a black hole. We will assume that the incident wave 
is monochromatic, long-lasting, and sufficiently weak that the gravitational field equations 
may be linearized. The scenario is then characterized by just two dimensionless parameters, 



(with units G — c — 1, used throughout). The coupling M\uj\ expresses the ratio of 
Schwarzschild horizon size r s to incident wavelength A, and < a* < 1 is a measure of 
the rotation rate of hole. In this paper, we concern ourselves only with the long- wavelength 
regime, in which M\u\ <C 1. 

The frequency u may take either sign, depending on the helicity of the incident wave. 
Positive io > corresponds to a circularly-polarized incident wave co-rotating with the black 
hole, whereas negative uj < corresponds to a counter-rotating helicity. If the scattering 
interaction is able to distinguish between these cases, then a polarization will result. 

Over the years, various authors [U EJ [H [7J [91 [TOj [EE] have shown that, for a non-rotating 
massive body (a* = 0) in the long wavelength limit [Mu <C 1) the cross section depends on 



M\u\ = nr s /\ 



and 



a, = a /M = J/M 2 
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the spin s of the scattered field as follows: 



. ttow\ s = 0, Scalar wave lal 

sm (0/2) L J 

s 2 ' 

sin T (6»/2) J ~~ 2' 

cos 4 (0/2) 



co ^ e ^ s=± Neutrino [b] 

' (2) 



1—77 1. Photon [c] 

snruTHT) s = 2, Grav. wave [d]. 



sin 4 (0/2) 

cos 8 (0/2) + sin 8 (0/2) 



The gravitational result is somewhat anomalous, in that it doesn't follow the same general 
rule [da/oin = M 2 cos 4s (6>/2)/ sin 4 (0/2)] as the other fields. This is related to the fact that 
the helicity of a gravitational wave is not conserved by the scattering process [8j [7J [11] . 

The influence of rotation on the scattering of low-frequency waves (M|o;| <C 1) has been 
considered by a number of authors p2J EE dH EEl EE E3 13 EJ [18]. It seems that no clear 
consensus has yet emerged. Considering electromagnetic radiation, Mashhoon [12] noted that 
"one must expect partial polarization in the scattered light when an unpolarized wave is 
incident on a Kerr black hole", and went on to observe that "the polarizing property of a 
Kerr black hole is probably maintained for very low frequencies" . However, the corollary that 
incoming gravitational waves will also be polarized at low frequencies was thrown into doubt 
by a study [7J in which the scattering amplitude was computed using Feynman-diagram 
techniques. The authors concluded that "the angular momentum of the scatterer has no 
polarizing effect on incident, unpolarized gravitational waves" (even though they found that 
unpolarized electromagnetic waves were polarized). Some doubts about the gauge- invariance 
of the results in [7] have been raised. To clarify the issue, the Feynman-diagram approach 
was recently revisited and improved in [18] . In a series of papers, Guadagnini and Barbieri 
[T5| [T6l [T71 [T8] have argued that the net polarisation V induced by rotating classical matter 
is 

^>0)-^<0) k 

% (oj > 0) + % (oj < 0) 
where 9 is the scattering angle, which is assumed to be small. 

In this paper we combine the partial wave approach of Matzner et al. [T9} [HI [2Ql [9] with 
the low-frequency asymptotics of Mano, Suzuki and Takasugi [23] to derive the lowest-order 
correction to cross section (E])[d] for a gravitational wave impinging along the rotation axis 
of a Kerr hole. We find a term proportional to au that couples the helicity of the incident 
wave with the spin of the hole. The total cross section is given in Eq. ( J72l) . The resulting 
polarization, given in Eq. (1731) . is in qualitative agreement with (J3j) at small angles, but differs 
by a factor of two. 

The remainder of this paper is organised as follows. In section 12.11 we briefly recap the 
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results of Matzner et al. [121 El EDI E] to write the cross section in terms of amplitudes which 
are expressed as partial wave series. In \2.2\ we define the phase shifts and discuss their 
asymptotic values [231 123] m the long-wavelength regime (M|a>| <C 1). In 12.31 we expand the 
relevant spin- weighted spheroidal harmonics to second order in au. In 12.41 we identify that 
part of the scattering amplitudes that is proportional to au, and compute the corresponding 
contributions to the cross section. We conclude with a brief discussion in section [31 



2. Analysis 

2.1. Partial Wave Series 

In the late 1970s, Matzner and co-workers [T9l El [201 EE] showed that the differential cross 
section for the scattering of gravitational waves by a rotating black hole can be written as 
the sum of the square magnitude of two amplitudes, 

% = w)?+\m\*. (4) 

where 9 is the scattering angle. In the special case in which the incident wave vector is parallel 
to the rotation axis, the amplitudes may be expressed as the following partial-wave series, 



_ oo 

9(fi) = - E E P (- 1 )' bM^L) - 1] -2Sf(0 ]a u) _ 2 S?(7r-e;au). (6) 



P=±l 1=2 

In these expressions, exp(2i5^ nu ) are phase factors to be determined from a radial equation, 
s 5 , ™(6 1 ; au) are spin-weighted spheroidal harmonics. Note the presence of the sum over even 
and odd parities, P — ±1. 

The first amplitude f(x) corresponds to (that part of) the interaction which preserves 
the helicity (i.e. for which the helicity of the scattered wave is the same as the helicity of 
the incident wave). The second amplitude g(x) corresponds to (that part of) the interaction 
which reverses the incident helicity. The helicity-reversing amplitude is non-zero because the 
phase shifts 5[ muJ depend on parity P. In this respect, gravitational wave scattering is unlike 
scalar, neutrino, or electromagnetic scattering. 

In [11] it was shown that, to lowest order in Mu, the Schwarzschild (a = 0) amplitudes 
may be written as 

f (x) - M C '* r(1 -' e) [|(1 + x)]2 m 

/ Schw (*) r(l + ie) -*)]!-* (7) 
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and 

^schw(x) = - x)] (8) 

where $ = 2eln |2e|, x = cos 8, e = 2Mu and T(z) is the Gamma function. In the following 
sections, we show that black hole rotation (a > 0) introduces a term in each amplitude which 
is proportional to au. 

2.2. Phase Shifts 

The phase shifts are found from solutions s Ri mLL) (r) to Teukolsky's [22] radial equation that 
satisfy the following boundary conditions, 

f £(trans) A _ s ifcr . 

s-n-Jmwl/V ~ S ( inc) - _. D (rcfl) _f2«+11 +iu;r ' ^ ' 

1 #L,, r L e mr * + JBL ,., V i^+xjg+wr, r* -> +oo 



Imuj ' Imu) 



Here, A = r 2 — 2Mr+a 2 , k = uo — a*m/ [2(1 + y/1 — a* 2 )], r* is the standard tortoise coordinate, 
s = —2 for a gravitational wave and -B^j, B^f/ and B^^ are complex constants. The 
phase shifts are found from 

^ Irau) 

where C is the Starobinskii constant, 

[Re(C)] 2 = ((A + 2) 2 + 4arwj - 4a 2 cj 2 ) [A 2 + 36amu - 36aV] 

+ (2A + 3)(96aV - 48aum) - lUu 2 a 2 , (11) 

and A = A\ m — 2mau + a 2 u 2 is the angular separation constant. 

Low-frequency analytic results for the phase shifts may be obtained via the formalism 
developed by Mano, Suzuki and Takasugi [231 [23] (MST). In the MST approach, reviewed 
in [25], solutions to the Teukolsky radial equation satisfying boundary conditions (Q are 
expressed as infinite series of special functions. Two series are used. The 'horizon' series of 2F1 
hypergeometric functions is convergent at all radii up to (but not including) spatial infinity. 
The 'far-field' series of Coulomb wavefunctions is convergent up to (but not including) the 
outer horizon. The complex constants Bj££j, Bj^j and B^^ are determined by matching 
the two series. It is found that 

,(inc) 



Dime, _ a ten 



sin yK[y + s — ie)) 



A+, (12) 

am ^vi if t o — «cn 

= Ae-cu- 1 " 2 * [K u {s) + i^K^s)} A v _. (13) 
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Here A s is just a normalization constant, and we follow the conventions of [23J by defining 



e = 2Muj and k = y/l - a* 2 . (14) 

Note that e must be positive in the MST expressions; the uo < results may be found 
via the symmetry of the radial function s Ri mu} = s R*- m -u from which it follows that 

T^(inc/refi) ^(inc/refl)* 

lm—u) l—muj 

The parameter v is known as the "renormalized angular momentum" and has the low- 
frequency expansion 



v = I + 



e 2 



21 + 1 



o * 2 , [(i + i) a -*T (i 2 -* a ) a 



+ 0(e 3 ).(15) 



1(1 + 1) (21 + l)(2Z + 3)(2/ + 3) (2Z-1)(2Z)(2Z + 1) 

Note the absence of a linear term in e. 

The coefficients K v and if_„_i may be computed via a complicated series expansion, 
detailed in [23]. This is not necessary for our purposes, because K_ v _i/ K u ~ 0(e 21 ^ 1 ). Hence 

,(inc) 
Imuj 



K_ p _i may be neglected, and the K v terms will cancel upon taking the ratio of Bt*£j and 



Al = r^-V^^-)e- ^f + | 5 + i€ \ £ <( S ), (16) 



n=—oo 



B (refl) 

The coefficients and A v _ are given by 

r(z/ + l + s-ie 

oo / i • \ 

A* = 2 - 1 - s e +-e- i ^^ +1+s ) e -^ 2 V (-l) n ) "TX OO, (17) 

n=— oo v 

where ~ (9(e' n '). Here, we wish to conduct an expansion accurate to second order in Mcv; 
hence we require explicit formulae for a u _ 2) ■ ■ ■ , a%. These were computed in [23] and are listed 
in 



Appendix A 



After inserting equations flA.lHA.5]) into (flOl) . and taking some care with series expansions 
in the small parameter e, we find 



exp(2z5r m J = e ^in2 ee -^ e -^-o gj + 1 ™) 4UW+D x + a ^ + Q(e 



,( 2 ) .2 , ^/ 3>, 



T(l + 1 + ie) 

for positive u. Note that there is no first-order term in the square brackets; the linear terms 
cancel exactly. The second-order coefficient is 
( 2 ) _ ima* 12ima* 
a ira ~ + ~ (Z-l)Z 2 (Z + l) 2 (Z + 2)' { } 

where m is the azimuthal number. Note that the e~ m ( v - 1 ) and e 4le / l ( l+l ) factors also give an 
/-dependent contribution at second order in e 2 . 
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Via the symmetry b[ v ^ c) = B^^ c) * it follows that exp(2iSf m _J = exp(2i5 l i mw )* . 
Hence the phase shift may be written more generally as 



exp(2^ m J = e ^\M e -^ e ~^-^/H ^ + 1 ie ) e ^/Ki+i) l + a Yy + (e A ) . (20) 

This expression is valid for either sign of e = 2Mu. 

In a previous study [11] , the low- Ma; approximations of Poisson and Sasaki [26] were 
used to show that the Schwarzschild (a* = 0) phase shifts are 

ex P (2^ m J = e 2 - ln2 -e- ^[| + |~g e^^) [l + 0(e 2 )} . (21) 

This is consistent with ( 1181) since re = 1 in the non-rotating case. Moreover, it is remarkable 
that, when a* = 0, equation ( |2TI) actually holds to one order higher in e than first supposed. 

2.3. Spheroidal Harmonics 

The spin-weighted spheroidal harmonics s S] n {9\ z = auj) are solutions of the equation 

1 d ( n dS\ ( 9 o „ m 2 2mscos9 „ 9 9 „ , \ „, „. 

sin 9— + [z 2 cos 2 9 - — y- r - r - 2zs cos - s 2 cot 2 + s + A lm }S(32) 



„( 2 ),2 



sin 9 d9 V ^ / V sin 2 sin 2 9 

where the spheroidicity parameter is 

z = aio. (23) 

These functions are regular at 9 = and 9 = ir, and are normalized so that 

J d(cos9) s Sr(9;z) s SJ l (9;z) = (24) 

Note that we suppress the azimuthal factor e im ^ throughout. 

The spheroidal harmonics may be decomposed into a sum over spherical harmonics 
s YJ n (9) of the same spin s and azimuthal number m. That is, 

oo 

a ST(0;z)= $J7V)> (25) 

j'=max(|m|,|s|) 

The expansion coefficients bj satisfy a five-term recurrence relation, 

(* 2 4 2 L) 4'4 + - 2-41.) 44 + (* 2 4S - 2-42 - *(* + d) 4° 

+ (A« +1 - 2z S c£> +1 ) 44 + ( A S +2 ) 44 = -*W-2> (26) 
where E im = A im + s(s + 1) and 



rffi s V fc m W cos .^(fl) = y (j, 1, m, 0|fc, m) (j, 1, -s, 0|fc, -s) , (27) 



4 2) = | ^ s y-(e)cos 2 e s rf(0) = i% + ^^i (j,2,m,0|fc,m) (j, 2, -s, 0|fc, -a) . (28) 
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The numbers (ji, j'2, mi, m^lj, m) are Clebsch-Gordan coefficients. For more details see, for 
example, Appendix A in [27]. 

We wish to expand the spheroidal harmonics to second order in z. Noting that 
&|2n ~ 0(z n ), ^ us m ake the expansion 

(29) 



6 (0 - z 2 d {0) 

u l~2 ~ * u -2' 

b ( {l 1 = zd®+z 2 d% 

bf =1+Z^\ 



6gi = + 

°Z+2 ~~ Z U +2- 



(0) 



.2^(1) 
+1' 



The normalisation condition implies that (do) 



unknowns {d_l, d_[, d_ 1: d Q 2 \ d+\, d+\, d+l} are determined from the equations 



(d®) 2 + (O 



(Oh 2 



(30) 
(31) 
(32) 
(33) 

. The remaining six 



(o) 
+1 



(l + l)d 
Id® 
2(1 + 



-sc 



(i) 



sc. 



(1) 



ZA W+1,M-1 + ^Iml^+l ^ c l+l, 



.(2) 



+ 1 



2(2/ + 3)dS 



.(2) 



J2) 
-'1+2,2 



2sc 



(i) 



(0) 



Z+2,Z+1 U +1) 



2(2/ l)di 2 — c l-2,l ^ SC J-2,J-1^-1> 

where the angular eigenvalue has been expanded in powers of z = ac<j as 



E lm = l(l + 1) + J2 E , 



Zm * 



(34) 
(35) 
(36) 
(37) 
(38) 
(39) 

(40) 



fc=i 



and the first few coefficients are determined by the identities 



E 



(i) 



lm 2SC U 



E {2) + c 



(1) 



(41) 
(42) 
(43) 



= 

(2) -2sc (1) d (0) -2sc (1) d {0) -0 
zz ZiC z,z-i"-i z ' b h,i+i a +i ~ u 

r(3) , (2) ,(0) (2) ,0 _ 9 (1) Al) _ 9 (1) ,(1) _ q 

For this calculation, we need only the m = 2, s = —2 harmonics. The required 
Clebsch-Gordan coefficients are listed in Appendix A , flA.6HA.13l) . Substituting the Clebsch- 
Gordon coefficients into ( 134H391 yields explicit expressions for the expansion coefficients 
{d {0 l d ( X d ( X d [ 2 \d% d% dfl}, which are listed in |Appendix Aj 

To compute the spheroidal harmonics explicitly we require expressions for the spin- 
weighted spherical harmonics _2^ 2 (a;), where x = cos 9. These may be found by acting 
on spherical harmonics of spin- weight zero, qY^x) =~ ' 21+1 



111 



Pi(x), with ladder operators [28] 
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The spin- weight is lowered with the operator 8, and the azimuthal number is raised with L + . 
These operators are defined by 

5 s Yr(x) = (VT^d x - ^=|) s Y l m {x) = -^(l + s)(l-s + l) s ^Yr(x), (44) 
L + s Yr(x) = - (VT^d x + S + mX ) s Yr(x) = v /(/- m )(/ + m + l) s ^+ 1 (x). (45) 



Here, d x is shorthand for the partial derivative with respect to x — cos 8. By acting with 
SL + SL + on o^°(x), it is straightforward to show that the spin-weighted harmonics in ([5]) and 
flU can be written 



W = V 4. (/-!)/(/ + !)(/ + 2) ' (46) 
t x P x {x) = (1 + x) 2 d x (l - x)d x d x (l - x)d x P l (x). (47) 

Their values in the forward and backward directions are particularly simple, 



*Y?(x = 1) = , _ 2 lf (a; = -1) = 0. (48) 

This implies that the values of the spheroidal harmonics in the forward direction are 



iST(x = 1) = J^+± (i + Sl z + S 2 z 2 + 0(z 3 )) (49) 



where 



Sl ~ JiTW (50) 

-3375 16 192(/ 2 + / + l) 2 32 

2 ~ (2/ + 3) 2 (2/- l) 2 + P{1 + l) 2 + (/ + l) 4 / 4 (Z + l) 4 / 4 ' ( ^ 

2.^. Scattering Amplitudes 

2. 4-1- Helicity- conserving amplitude Let us first consider the helicity-conserving amplitude 
/ defined in (JHJ). To begin, we note that the sum of the positive and negative-parity phase 
terms can be written 

I + e ™LJ\ = e ^x r(/-l~*e)r(Z + 3) r (2) 2 3 i 

2v e +e 7 6 r(/ + 3 + ze)r(/-i) I +Ule j J ' 

where e %x = e 2leln \ 2e \ e - teK anc i the second-order coefficient is 

flS^^fe^-^ + TTT^- T^ , ,1 » ■ (53) 



ioji v e 2 



'/(Z + l) /(/ + 1) (2/ + 3)(2/-l)' 
Inserting the expansions of the spheroidal harmonics (125} l4~9l) , and using fl52|) , we may write 
the / amplitude as 



e 



= 2^ F(X) ' (54) 
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where 

F(x) = |>Z + l) r ( ( j + 3 + £) (l + AS^ 2 + ^(e 3 )) (1 + + + 0(, 3 )) ^, (55) 

and 



" 2/ + l(/-3)(/-2) ^ 2 v ; V 2Z + l(Z-2) z_1 u ; ' v ; 
2Z + 3(Z-l) t(0 D , . /2Z + 5 (Z-l)Z , {0 



+ ^a + i^^ 1 ^ + V2TTI (/ + 3)(z + 4) ^^ ^ 

and are the expansion coefficients defined in fl25l) . 

To compute the higher-order terms in the sum F(x), we substitute in the explicit forms 
for the coefficients calculated in Section |2T3"1 Next, we rewrite the sum so all terms have a 
common factor of Pi. To demonstrate this process, let us begin by considering just the linear 
term in u, which we denote F(^. We find 




f m (x) = a^}_^(2i + | ^ ~ ^ \ || + ^ ifflfu*) + s^x) 



oo 
1=2 

,2 




= o + 

It is straightforward to verify that the terms in parantheses on the second line cancel exactly, 
and so a linear-in-u; term in F(x) is not present. 

Note that to obtain Pi(x) as a common factor in (1571) we redefined the summation variable 
(i.e. Z — > Z ± 1). Care must be taken in this process, since it changes the lower limit of the 
sum. In the above, we added an 'extra' P2(x) term to the series without consequence, because 
the coefficient d^\(l — 1) is zero for I = 2. We also neglected a P\(x) term (present in the 
top line but not subsequently). This is justified for our purposes because, to obtain the final 
amplitude /, we act on F with C x , and C x P\{x) = 0. 

Let us now repeat this process and focus only on the part which is linear in a. This time, 
we will keep terms up to second order in u. Such terms arise from (i) the a*-dependent part 
of the phase /^e 2 , and (ii) the linear term S\ coupled to an e term arising from the effect 
of redefining the summation variable. We may split the result into two parts: a sum which 
turns out to be zero, and an Z = 2 term which arises from the change of summation variable. 
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That is, 

F(x) « F (x) + ze (f^ + if =2) ) + C(^ 2 , e 2 ) + 0(u; 3 ). (58) 
The lowest-order term 

ft = E(2« + l)g^|flW (59) 
z=2 v 1 1 ; 

was defined and summed in [11]. The polarizing terms are 

r(Z - 1 - ze) 



F 1 (E) (x)=^^(2/ + l) 



IV • 3 • «) 



and 



-2(/ + l) (21 + 3) (0) 4 2/ {2l-\) {0) \ 

(/-!)(/ + 3) V(^Tl) - l( + } " + + (/-2)(/ + 2)V(2TTl) rf+l(/ ~ 1} J P ' (X) 
(60) 

^ 2) (x) = -«^gft(x). (61) 

(7=2*) 

The term zei^ (x) is responsible for the lowest-order polarizing effect. It leads to a term 
in the scattering amplitude which depends on the sign of to. That is, / ~ /schw + / po i where 



/poi = -2iWW* 



;i + O(e)) • (62) 



r(l + ie) 

Since / po i is in phase with /g c h w , this gives a first-order contribution to the cross section, 



r2 cos 8 (0/2) 



2 |/schw/poi| = ~ 4 ^ M — ITg^y - (63) 

2.4-2. Helicity-reversing amplitude We now repeat the analysis for the helicity- reversing 
amplitude, g(9). To simplify matters, we only expand the relevant terms to first-order to 
recover the polarizing correction. Let us begin by writing the amplitude as 



g(9) = —£- x G(x), wh ere G(x) = ]T(-l) z WW(i + S lZ + 0(z 2 )), (64) 

1=2 

and 



T(l -I - it) V(l - 1) / 2a*m ., 



Wl = 6 - r(z + 3 + ,e[r(z + 3) / { 1 + WTTf + °^ ) • ™ 



^ = fl(- x ) + ^ ____ 2 dW(0^-i(-x) + ^ 2/ - n - 3 < u l(/)P m (-x) + 0(z 2 ). (66) 

As before, we may move terms of the series up or down (I — > i± 1) to obtain a common factor 
of Pi(—x). As before, we find that, to lowest-order, the terms in the sum cancel out, leaving 
only an I = 2 term which arises from redefining the summation variable. That is, 

G{x) « G (x) + ze (cf } + G { l =2) ) + 0(z 2 , e 2 ) + O(co 3 ), (67) 
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and 

G™(x) - 6:fv IV ^- 1 -") r C - V ((21 + 1)^ 2 ( / + 2 ) 2 + 2 ( f - 1 H 

Gl w " b *2^ ^ r(/ + 3 + *e) r(/ + 3) V V + 1 ) /2 (^ n2 ' 



«=2 



= 0, (68) 

Gf =2) W=2 ? ^^4p 2 (-x). (69) 
1 (5 + zej 

The zeG\ (x) term is responsible for the lowest-order polarizing effect. It leads to a term 
in the scattering amplitude which depends on the sign of uj. That is, g ~ gschw. + 9 P oi, where 



g po \ = 2Maue lx 



-(1 -x) 
2 V ' 



1 + 0(e)) . (70) 



r(l + ze 

The first-order contribution to the helicity-reversal cross section is 

2\g* Schw9pol \=4aujM 2 sm 6 (e/2). (71) 

3. Discussion and Conclusion 

In the preceding sections we have shown that, in the low- frequency regime (M|u;| < 1), the 
scattering cross section for a gravitational wave impinging along the axis of a rotating black 
hole is 

- Si i 1 - - *w>] + 4 - ^ ■ (72) 



The first term in (172]) arises from the helicity-preserving interaction; the second term arises 
from that part of the interaction that reverses the helicity of the incident wave. 

The cross section (I72p clearly depends on the sign of u, so scattering from a rotating black 
hole induces a partial polarization in an initially unpolarized beam. The total polarization is 

S(">0)-S("<0) = _ 4aU sia , (m) <W(»/2)-»n'(»/2n 



" - gfc > oi + jfc < ,i ^ - 4 ^' sta '<^> U^iw^S J • < 73 > 

In the small-angle limit, the polarization is exactly half that predicted by Guadagnini and 
Barbieri [151 ESJ El EE] for scattering from classical rotating matter (eq. [3]). 

There are a number of possible extensions of this work. Firstly, the analysis could be 
extended to encompass waves approaching at arbitrary angles of incidence with respect to the 
rotation axis. The relevant partial wave formulae, given in [9], involve an additional sum over 
the azimuthal numbers m. We would expect to find maximal polarization for waves impinging 
along the axis of rotation, and zero polarization for waves approaching in the equatorial plane. 
Secondly, the analysis could be repeated to compute the polarization of electromagnetic waves. 
Whilst we believe that the main features of this analysis follow through without difficulty, it 
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remains to be checked. Finally, in a complementary work [29], numerical methods have been 
employed to compute cross sections for arbitrary wavelengths. Excellent agreement between 
numerical and analytical results is found in the low-frequency regime. In addition, interesting 
higher-order effects are observed for general Mu, such as interference fringes and glory halos. 
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Appendix A. Coefficients 



The coefficients a v n appearing in the MST formalism ffTTl - fT6]) are 



(/ — 1 + s) 2 (l + s) 2 [(I — 1)k — ima*} [Ik — ima,] 2 



4(1 - l)/ 2 (2/ - 1) 2 (2/ + 1) 



t l + 0(e 



.(/ + s) 2 [Ik — ima*} (I + s)' 



2/ 2 (2/ + r 



2/2(2/ + 1) 



l + i- 



(l + l-s) 2 [(/ + 1)k + ima*] 
2(1 + 1)2(2/ + 1) 



e + 



Ik — ima* 
(/-l)/ 2 (/ + l) 

(l + l-s) 2 



ma*s 



(A.l) 



e 2 + 0(e 3 ), (A.2) 



2(1 + 1) 2 (2/ + 1) 



.(/ + 1)k + ima* 2 
1 — -7—— —ma*s 



/(/ + l) 2 (/ + 2) 



1/ 

«2 



(/ + 1 - s) 2 (l + 2 - s) 2 [(I + 1)k + ima*] [(/ + 2)k + ima 



*\ ,2 



e 1 + 0(e" 



(A.3) 
(A.5) 



4(/ + l)2(/ + 2)(2/ + l)(2/ + 3) 2 
With the Cordon-Shortley phase convention, the Clebsch-Gordan coefficients appearing 
in ([34H39D are 



(/, 1, 2, 0|/ + 1, 2) = [(/ -l)(l + 3)/(2I + 1)(/ + 1)] 1/2 

(/,l,2,0|/,2) =[4//(/ + l)] 1/2 , 

(/, 1, 2, 0|/ - 1, 2) = - [(/ - 2)(Z + 2)/(2l + 1)/] 1/2 , 



and 



(/,2,2,0|/ + 2,2) 

(/,2,2,0|/ + 1,2) 
(/,2,2,0|/,2) 



(2/)! (/ + 4)! (/! i 1/2 
_"(2/ + 4)!(/ + 2)!(/-2)! 

12 (/-!)(/ + 3) 1 1/2 
(2/ + l)/(/ + l)(/ + 2)J ' 
-2(/-3)(/ + 4) 
[(2/- l)(2/)(2/ + 2)(2/ + 3)] 1/2 ' 



(A.6) 
(A.7) 
(A.8) 

(A.9) 

(A.10) 
(A.ll) 
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<Z,2,2,0|Z - 1,2) =- 



1/2 



(Z,2,2,0|Z-2,2> 



12 (Z-2)(Z + 2)" 
(2Z + 1) (Z-l)/(/ + l)_ 

(2Z - 3)! (/ + 2)! (Z - 2)P 1/2 



(2Z + 1)! Z! (Z-4)! 
The expansion coefficients {d_l, d^\, d^\, d!^\ d+\, d+\, d+l} appearing in ( J3~3l are 



d% 



7(0), 



2[(2Z + l)(2Z + 3)]-^ ^-^t^ . 



-2[(2Z + 1)(2Z- 1)] 



_ 1/2 (Z-2)(Z + 2) 



4[(2Z + l)(2Z + 3)] 



_ 1/2 (/_!)(/ + 3)((Z + l) 2 - 8) 



Z(Z + l) 4 (Z + 2) 

72 



I f/9I , 1V9I oM-i/2 (Z-3)(Z-2)(Z + l)(Z + 2)(Z-8) 
-L^ + lj^ 3)j (2Z - 1)2(Z - 



(A.12) 
(A.13) 

(A. 14) 

(A.15) 
(A.16) 

(A.17) 

(A.18) 

(A.19) 
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